Incremental Redundancy Cooperative Coding for Wireless Networks: 
Cooperative Diversity, Coding, and Transmission Energy Gain 

Ruoheng Liu, Predrag Spasojevic, and Emina Soljanin 



' Abstract 

o ■ 

I We study an incremental redundancy (IR) cooperative coding scheme for wireless networks. To exploit the 

, distributed spatial diversity benefit we propose a cluster-based collaborating strategy for a quasi-static Rayleigh 

i fading channel model. Our scheme allows for enhancing the reliability performance of a direct communication 

■ over a single hop. The collaborative cluster consists of M — 1 nodes between the sender and the destination. The 
transmitted message is encoded using a mother code which is partitioned into M blocks each assigned to one of 
M transmission slots. In the first slot, the sender broadcasts its information by transmitting the first block, and its 
helpers attempt to relay this message. In the remaining slots, each of left-over M — 1 blocks is sent either through 

O , ■ a helper which has successfully decoded the message or directly by the sender where a dynamic schedule is based 

on the ACK-based feedback from the cluster By employing powerful good codes including turbo, low-density 
] parity-check, and repeat-accumulate codes, our approach illustrates the benefit of collaboration through not only 

OO 

l/^ . a cooperation diversity gain but also a coding advantage. The basis of our error rate performance analysis is the 

o 

■ threshold behavior of good codes. We derive a new simple code threshold for the Bhattacharyya distance based 



o 
o 



on the modified Shulman-Feder bound and the relationship between the Bhattacharyya parameter and the channel 
capacity for an arbitrary binary-input symmetric-output memoryless channel. The study of the diversity and the 



C/2 

Q . coding gain of the proposed scheme is based on this new threshold. An average frame-error rate (FER) upper bound 

and its asymptotic (in SNR) version are derived as a function of the average fading channel SNRs and the code 
' threshold. Based on the asymptotic bound, we investigate both the diversity, the coding, and the transmission energy 

CZ ^ gain in the high and moderate SNR regimes for three different scenarios: transmitter clustering, receiver clustering, 

and cluster hopping. Furthermore, given a geometric distance profile of the network, we have observed that the 
energy saving of the IR cooperative coding scheme is universal for all good code families. 
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I. Introduction 

To overcome fading, wireless networks employ various diversity techniques, e.g., channel interleavers, multiple 
antennas, frequency hopping, etc. In [1], Sendonaris, Erkip, and Aazhang have proposed the so-called user- 
cooperation diversity where users partner in sharing their antennas and other resources to create a virtual array 
through distributed transmission and signal processing. Only limited time and frequency diversity is available when 
the transmitted signal experiences a quasi-static frequency-flat fading and the receiver has a strict decoding delay 
constraint (as in speech and video transmission). In this case, introducing spatial diversity through cooperation is 
especially beneficial. 

Motivated by the increasing importance of cooperative communication in wireless networks, a large of valuable 
work has gone into designing the collaborative protocols, assessing the information-theoretic limits, and evaluating 
the cooperative benefit in the past few years. An information-theoretic analysis of several cooperative protocols has 
been reported in [2], [3] based on repetition-coding and space-time-coding. In an independent, relevant work [4], 
the authors have studied the outage probabihty and the asymptotic error probability performance calculations of a 
two-user incremental redundancy coding scheme based on Gaussian codebooks. An practical coding strategy for 
two-user collaborating transmission based on rate compatible punctured convolutional codes has been considered 
in [5], [6]. Cooperation schemes employing distributed turbo codes have been studied in [7]-[9], where [8] has 
focused on the asymptotic error performance analysis. 

In this paper, we study an incremental redundancy (IR) cooperative coding scheme for a multiple-helper 
wireless network. To exploit the spatial diversity benefit we propose a cluster-based collaborating strategy for 
a quasi-static Rayleigh fading channel model and based on a network geometric distance profile [10]. The cluster- 
based collaborating strategy enhances the reliability performance of a direct communication over a single hop. A 
collaborative cluster consists of M — 1 nodes assisting the sender in a two-hop strategy. The transmitted message 
is encoded using a mother code which is partitioned into M blocks each assigned to one of M transmission slots. 
In the first slot, the sender broadcasts own information by transmitting the first block, and its helpers attempt to 
decode the message. In the remaining slots, each of the left-over M — 1 blocks is sent either through a helper 
which has successfully decoded the message or directly by the sender based on an acknowledgement (ACK) driven a 
dynamic schedule. By employing the powerful good codes [11], [12], e.g., turbo, low-density parity-check (LDPC), 
and repeat-accumulate (RA) codes, whose performance is characterized by a threshold behavior, our approach 
illustrates the benefit of collaboration through not only a cooperation diversity but also a coding advantage. 

To investigate benefits of the proposed IR cooperative coding scheme, we evaluate its frame error rate (FER) 
performance for a quasi-static frequency-flat Rayleigh fading channel. The basis of our error rate performance 
analysis is the so-called code threshold of a good code ensemble. Our focus is on the code threshold for which 
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the maximum likelihood (ML) decoding word error probability vanishes whenever the channel quality is above 
the threshold. In fact, we propose a new Bhattacharyya distance simple code threshold for parallel channel 
transmissions based on the modified Shulman-Feder reliable channel region [13]-[15] and the relationship between 
the Bhattacharyya parameter and the channel capacity of an arbitrary binary-input symmetric-output memoryless 
(BISOM) channel [16]. The derived threshold characterizes a reliable communication condition for parallel channel 
transmissions with a single constraint in terms of the average Bhattacharyya parameter. Compared with the union 
Bhattacharyya (UB) code threshold introduced in [17], the new code threshold is almost tighter 1 dB for an example 
turbo code transmitted over a binary-input AWGN (BI-AWGN) channel, as illustrated in Example [T] 

Next, based on the simple code threshold and the outage concept, we derive a PER upper-bound which predicts 
well the fading channel simulation results for the proposed cooperative coding scheme. Compared with the result 
in [6] (and the block-fading channel PER bound in [18]), the derived PER bound is a function of the single code 
threshold parameter instead of a family of weight enumerators. Hence, this simple PER bound is not only easy 
to compute but is also insightful in that it allows for a large SNR asymptotic analysis of the cooperation scheme. 
In this paper, closed-form asymptotic PER upper-bounds are obtained for three different scenarios: transmitter 
clustering, receiver clustering, and cluster hopping. These bounds allow for illustrating cooperative diversity and 
coding gains in the high SNR regime. Pinally, we express PER bounds in terms of the distance profile and investigate 
cooperation benefits in terms of the energy efficiency gain for different collaborative cluster sizes and normalized 
cluster distances. 

The paper is organized as follows. We describe the system model in Sec. JI] We study the threshold behavior of 
good code ensembles over BISOM channels in SecJII] An IR cooperative coding scheme is introduced in Sec. IIVI 
We derive an upper bound on the scheme PER in Sec. |Vl and its asymptotic versions for different cooperation 
scenarios in Sec. |Vll Simulation results and the collaborative cluster design are discussed in Sec. IVIII Pinally, we 
summarize our results in Sec. IVIIII 

II. System Model 

In this section, we introduce a collaborative cluster and describe a fading channel model based on the geometric 
distance profile of the network. 

A. Collaborative Cluster 

As shown in Pig. [T] we consider a single sender-destination pair and several intermediate nodes. A group of 
M — 1 nodes forms a collaborative cluster S which assists the sender in a two-hop scheme. The sender serves as the 
regional broadcast node, and each cluster member, termed helper, attempts to relay the package to the destination. 
Por notational convenience, let Node denote the sender. Node m, for m = 1, . . . , M — 1, denote a helper, and 
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destination 



collaborating cluster 



Fig. 1. Cooperative routing in wireless network with a geometric distance profile 
Node M denote the destination. 

Remark 1 The proposed cluster-based collaborating strategy can be embedded into an existing noncooperative 
route. For example, the routing metric most commonly used in existing ad hoc routing protocols is the minimum hop- 
count, which results in a long distance traveled by each hop [19]. Hence, there may exist several intermediate nodes 
between the (hop) sender and the (hop) destination for each non-cooperation hop. This is one of the motivations 
for the IR cooperative coding scheme proposed in this paper, namely, the cooperative scheme is a supplement to 
current routing algorithms for overcoming the "dynamic" deep fade. 

B. Channel Model 

We consider the quasi-static frequency-flat Rayleigh fading channel model [20] where the fading coefficient is 
random but invariant during the transmission interval T. The discrete-time channel model is 

yi,j = ^hj + ^hj for i G {0, . . . , M - 1} and j G {1, . . . , M} (1) 

where is the signal transmitted by Node i, L is the path loss exponent, dij, aij, and Zij are the distance, 
fading coefficient, and background noise between Nodes i, j, respectively, and yij is the signal at Node j received 
from Node i. We will focus on the transmitted signal alphabet is binary, i.e., Xj G {+VE, —Ve}, where E is the 
transmitted symbol energy which is identical for all cluster nodes and the sender, and the symbol "+-\/E" (and 
"—■s/E") represents the coded symbol "0" (and "1"). We assume that Zij is modeled as the mutually independent 
additive Gaussian noise J\f{0, 1/2), and Uij = la^jp is the exponentially distributed "channel power" with mean 
1. Then, the average and instantaneous SNRs of the signal at Node j received from Node i can respectively be 
expressed as 

SNRij^E-d-f, and 9i,j ^ i^ij ■ SNRij . 
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We further assume that decoding is done with the knowledge of the fading coefficients. Note that, for a given fading 
coefficient a^ j, the channel ([T]) is a BI-AWGN channel. This motivates us to study the error performance of binary 
good code ensembles transmitted over BISOM channels in the next section. 



III. Threshold Based Performance Analysis for BISOM channels 

In this section, we introduce the notation and the preliminary material on BISOM channel measures and the 
weight spectrum of code ensembles. Next, we study the threshold behavior of good code ensembles and derive a 
new simple code threshold for parallel channel transmissions. Results given in this section are used as the foundation 
for the error performance analysis of the IR cooperative coding scheme. 



A. Channel Capacity and Bhattacharyya Parameter of BISOM Channels 

Here, we consider a binary input memoryless channel with the output alphabet y and transition probabilities 
p(y|0) and p{y\l), y £ y. We say that the channel is symmetric if p{y\0) = p{—y\l). We first study the channel 
capacity and the Bhattacharyya parameter for BISOM channels. The Bhattacharyya parameter is widely used to 
characterize the "noisiness" of the channel (e.g., see [15], [17]). The capacity of a BISOM channel is achieved by 
the uniform input distribution and can be expressed in terms of p{y\0) as follows 



c = i 

2 



■yey 



y&y 



p{y) 



log; 



p(y|o) 







p{Y) 

p(y|o)+p(-y|o) 



log 



p{Y\Q) 



Similarly, the Bhattacharyya parameter 7 is 



(2) 



7 = 5Z VP(y|Orai) 

y&y 

We also consider the cutoff rate for the BISOM channel. 



p{-Y\0) 



p{Y\0) 



(3) 
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l-log2(l+7)- 



(4) 



Now, we establish a general relationship among these three information-theoretic quantities in the following lemma. 



Lemma 1 Let B = l — j denote the Bhattacharyya rate. For a BISOM channel with a transition probability p{y\0), 
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the channel capacity, the Bhattacharyya rate, and the cutoff rate satisfy 



C>B>Rq. 



(5) 



Proof: The proof in Appendix O follows the Jensen's inequality. ■ 

In Lemma [T] we propose a new channel quality measure called Bhattacharyya rate which is between the channel 
capacity and the cutoff rate. In particular, the Bhattacharyya rate is equal to the channel capacity for the case of 



- ■ Capacity 

— Bhattacharyya rate 
- Cutoff rate 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 
erasure probability 

a. BEC 




-15 -10 -5 

SNR (dB) 

b. BI-AWGN channel 



Fig. 2. Capacity, Bhattacharyya rate, and cutoff rate for BECs and BI-AWGN channels 



the binary erasure channel (BEC). For the case of a BI-AWGN channel with the SNR A, 



1 f+oo r- I— 

C{\) = 1-^1 e"(s^-^) log2 (1 + e-^y^) dy 



B{X) = l-e"^ 
i?o(A) = l-log2(l + e-^) 



(6) 

(V) 
(8) 



Clearly, for a BI-AWGN channel, the numerical calculation of the Bhattacharyya rate is much easier than the 
computation of the channel capacity which requires integration form from —oo to oo. The channel capacity, the 
Bhattacharyya rate, and the cutoff rate are illustrated in Fig. |2] for BECs and BI-AWGN channels. 



B. Weight Spectrum Properties of Good Code Ensembles 

In this paper, we consider good binary linear code ensembles [11], [12] whose performance is characterized by 
a threshold behavior. The class of good codes includes turbo, LDPC, and RA codes, and some recent variations. 
Due to the outstanding performance of such random-like codes, we employ these in our IR cooperative coding 
scheme. In order to evaluate the error performance of the cooperation scheme, we study code thresholds based 
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on the weight spectrum of code ensembles. The weight spectrum for various good code ensembles were studied 
in, e.g., [17], [21]-[23]. In this subsection, we briefly introduce some notation on weight spectra of binary linear 
code ensembles, review the UB code threshold proposed in [17], and describe a relationship between the UB code 
threshold and the code rate. 

For a given binary linear code ensemble [C] of length A'^ and rate R, let Af^^ denote the average number 
of codewords of Hamming weight h, termed the average weight enumerator (AWE). Let 6 = h/N denote the 
normaUzed Hamming weights, then the asymptotic normalized exponent of the weight spectrum with respect to the 
codeword length is defined as 

r['^](5) = limsupi^^4^ (9) 

where the superscript [C] denotes a binary code ensemblq^. To simplify notation, hereafter, we drop the superscript 
[C] from code parameters, e.g., in A^\5) and when we do not specify a particular code ensemble. Now, the 
Shulman-Feder (SF) distance [13] 

^ = sup [r{5) - r[^^]((5)] logs e (10) 
0<5<1 

measures the weight spectrum distance between the random binary code ensemble [JZB] and the code ensemble 
[C], where 

r'^'^^\5) = H{5) - {I - R)\n2 
is the asymptotic normalized exponent of the weight spectrum for the random binary code ensemble and 

H{5) = -6 ln{6) -{I -6) ln(l - 6) 
is the binary entropy Junction. Following [17], we define the UB code threshold of a code ensemble [C] as 

co^sup^. (11) 

0<<5<1 

In this paper, we consider a family of good code ensembles whose weight spectra satisfies the following condition: 
1) For a given binary code ensemble [C], there exists a sequence of integers D]y{< N) such that —>■ oo and 



D 



lim VA = 0; (12) 



h=l 

2) and the UB code threshold cq is finite. 

Finally, we establish a relationship between the code rate R and the threshold cq in the following lemma. 

'We shall use the symbol C to denote a binary codeword, {C} to denote a codebook, and [C] to denote a binary code ensemble. 
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Lemma 2 For a good binary linear code ensemble [C] of rate R, the UB code threshold cq is lower bounded by 
the following function of R, 

Co > -ln(l - i?). (13) 
Proof: We show the proof in Appendix |D] ■ 



C. Threshold Behavior of Good Codes for BISOM Channels 

The basis of our analysis is the threshold behavior of good codes. Here, we first review the work in [17], [24], 
which has studied the error performance of good code ensembles (e.g., turbo codes) based on the UB code threshold 
Cq. Next, we introduce a tighter Bhattacharyya distance code threshold (compared with cq) under ML decoding. 
Based on the new threshold c*, we derive a coding theorem for good code ensembles whose transmission takes 
place over independent parallel channels. 

1) UB Code Threshold: In [17], Jin and McEliece have shown that, for a binary-input memoryless channel, if 

7 < exp(-co), (14) 

the average ML decoding word error probability approaches zero. Inequality ([T4l ) describes a reliable communication 
condition for the code ensemble [C] based on a single channel parameter 7 and a single code parameter cq, where the 
Bhattacharyya parameter 7 represents the "noisiness" level of the channel and the UB code threshold cq characterizes 
weight spectrum properties of the code ensemble. This result is based on the classical union bound. Hence, the 
threshold cq is loose. Based on cq and the random assignment method (described below), we have studied this 
bound for parallel BISOM channels [15], [24]. 

2) Simple Code Threshold: Following the approach in [14], for a given good code ensemble [C] of rate R, we 
partition the normalized Hamming weights 6 into two disjoint subsets, 

^^J{P) ^ {6 : < S < 0.5 - P or 0.5 + P < 6 < 1} 
^'^(P) ^ {5 : 0.5 - P < (5 < 0.5 + P}, for P G [0, 0.5). 

Next, we define a new Bhattacharyya distance code threshold by optimizing the weight partition parameter P as 
follows 



= mill {cp : cp > — ln(l — R — 

0<P<0.5 



(15) 
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where 



CP A (16) 

<5e*(P) 



denotes the restriction UB code threshold corresponding to the weight subsets ^(P), and 

Cp=< 



sup [r{5) - r[^^l(5)] loga e, P > 
5e*<^{P) (17) 



0, P = 0. 

denotes the restriction SF distance corresponding to the weight subsets "^^{P). Note that, in the case P = 0, 

cp = Co and = 0. 

Hence, Lemma |2] and ([TSl l imply that the new code threshold c^, < cq. 

Now, we consider the average error rate performance of a code ensemble [C] transmitted over Q parallel channels. 
Following the random assignment approach [24], we assume that the bits of the transmitted codeword are randomly 
assigned to parallel channels so that each bit is transmitted over Channel j with the a-prior probability tj, where 
Ylf=o '^j — 1 ^^'^ '^j > for j = 0, . . . , Q — 1. We have the following parallel channel coding theorem for the 
code ensemble [C] based on c*. 

Theorem 1 Let us consider a linear binary code ensemble [C] whose transmission takes place over Q independent 
parallel BISOM channels. Assume that the bits of the transmitted codeword are randomly assigned to the channels 
with assignment rates {tq, . . . , tq_i}. Let 

Q-i 
j=0 

be the average Bhattacharyya parameter over Q parallel channels, where 7j is Bhattacharyya parameter of Channel 
i,fori = 0,...,Q-L If 

7 < exp(-c^), (19) 

then the average ML decoding word error probability Pwi^,N) ^ 0. 

Proof: The proof is based on the relationship between the channel capacity and the Bhattacharyya rate in 
Lemma [U the lower bound on the UB code threshold in Lemma |2j and the modified Shulman-Feder reliable channel 
region in Appendix lAl (please also refer to [12] for details). We provide the proof of Theorem[T]in AppendixlE] ■ 

In Theorem[T] the reliable communication condition ([T9l l is a simple constraint in terms of the average Bhattacharyya 
parameter. The new code threshold allows for characterizing more complex coding schemes (e.g., communication 
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over block-fading channel where the error performance requires averaging over all possible channel realizations) in 
a simple and more accurate manner. Hence, we refer to as the simple code threshold for [C]. The above simple 
code threshold theorem describes the asymptotic result where we let the codeword length N tend to infinity. The 
following example illustrates that the simple code threshold aids in estimating the error performance of long codes 
with fixed codeword lengtlj^ 



Example 1 (simple threshold for turbo codes) Here, we study UB and simple code thresholds of a R = 1/7 
turbo code. The turbo encoder consists of J = 3 recursive convolutional encoders with two random interleavers. 
The component code transfer functions are 
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Fig. 3. UB and simple code thresholds for an AWGN channel (turbo code of rate i? = 1/7 and length N — 5376, 10752) 



Gi = (1, 13/15, 17/15) and G2 = G3 = (13/15, 17/15). 

We compute the AWE based on the technique in [22]. By applying (fTTI ) and ( [T5] ), we calculate the UB threshold 

Cq'^''^ PS 0.21 and the simple threshold c^^'^^ « 0.17. As shown in Fig. [3l we compare UB and simple thresholds 

with simulation results for codeword length = 5376 and N = 10752 under iterative decoding when the turbo 

codes are transmitted over a BI-AWGN channel. Fig. |3] illustrates that the cliff of WER curves becomes sharp as 

the codeword length N increases. We observe that the simple threshold predicts well the cliff of the simulated 

word error probabiUty, and that the gap between clf^^ and ci'^'"' is almost 1 dB. We also consider the capacity Umit 

of the BI-AWGN channel. To this end, we set C(A) in Q equal to the code rate R, and determine the capacity 

"^For typical practical systems, the codeword length is fairly large [25], e.g., in CDMA2000 standard [26], the encoder allows for a variable 
input length up to A' ~ 20730. 
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threshold SNR. Fig. [3] illustrates that the capacity threshold may not predict well the error performance of turbo 
codes. 

3) Punctured Code Threshold: For punctured codes, we may assume that punctured bits are sent to a dummy 
memoryless component channel whose output is independent of the input, i.e., 7p = 1, whereas, non-punctured 
bits are transmitted over the real channel with the Bhattacharyya parameter 7. Let the puncturing rate be 1 — r, 
the average Bhattacharyya parameter is 7 = 7 • r + 1 • (1 — r). Hence, we can rewrite the reliable communication 
condition ( fT9l ) as 

r 

Analogous to Theorem [T] we have the following result for a (randomly) punctured code ensemble. 

Theorem 2 Consider a good codes ensemble [C] with a finite code threshold defined in (ITSl) . Assume that the 
coded symbols are randomly and independently punctured, so that each bit is punctured with a-priori probability 
(punctured rate) 1 — r. If r > 1 — exp(— c^,), there exists a punctured code threshold 

X(r) = In ^- (21) 

exp(-Cj,) - (1 - r) 

such that, if 

-ln7>x(r), (22) 

the average ML decoding word error probability approaches zero as ^ 00. 

Proof: The proof follows in a straightforward manner from Theorem [T] and (|20t . ■ 



We note that the left hand side of (|22i) is the channel Bhattacharyya distance, which is equal to the received SNR 
for a BI-AWGN channel. Thus, the simple code threshold is an SNR threshold for an AWGN channel, i.e., if the 
received SNR is larger than the punctured code threshold x{t)^ the average ML decoding word error probabiUty 
decays to zero as the codeword length N tends to infinity. We show the SNR threshold behavior of punctured turbo 
codes in the following example. 

Example 2 (threshold for punctured turbo codes) We consider an example punctured turbo code transmitted 
over a BI-AWGN channel. The mother turbo code of rate R = 1/7 and length N = 26880 is described in 
Example [T] We study the word error probability performance of its punctured versions of rate Rp = 1/5 (by 
setting T = 5/7) and Rp = 1/3 (by setting r = 3/7). Based on the simple threshold ci 0.17, we calculate the 
punctured code thresholds x(5/7) ~ 0.25 and x(3/7) ^ 0.45. As shown in Fig. |4l we compare the simulation results 
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Fig. 4. UB and simple code thresholds for an AWGN channel (mother turbo code of i? = 1/7 and N = 26880) 



with punctured thresholds and capacity SNR limits corresponding to C = 1/5 and C = 1/3. Fig. |4] illustrates that 
our punctured code thresholds are in a very good agreement with the numerical value observed in the simulation. 

Remark 2 Self-Decodable Condition: Theorem [2] illustrates that if 

r > 1 - exp(-c*), (23) 

then the punctured code threshold x{'^) exists for a good mother code ensemble [C]. It implies that ( [23l ) is a 
necessary condition which guarantees that a randomly chosen code (sequence) from the punctured code ensemble 
is self-decodable with probability one. Hence, we refer to (1231 ) as a self-decodable condition. 

IV. IR Cooperative Coding Scheme 

In this section, we introduce an IR coding collaboration scheme. We assume that the sender and helpers acquire 
a common time reference and operate in a single radio frequency band W . The medium access control (MAC) 
scheme is based on a time-division scheme. The transmission interval T is partitioned into M non-overlapping 
slots of duration tqT, . . . , tm-iT, where Tj, for j = 0, . . . , M — 1, is referred to as the assignment rate for slot j, 
X^jlo^ Tj = 1, and the time allocation {tq, . . . , tm-i} is predetermined for all nodes. 

Each node in the wireless network has an encoder (associated a codebook {C} of length N and rate R), a 
decoder, and a mapping device. The sender (Node 0) encodes the information and obtains a mother codeword C. 
The mapping device partitions the codeword C into M blocks corresponding to M non-overlapping slots. Block j 
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of length TjN is transmitted in slot j = 0, . . . , M — 1. For analysis simplicity, the partition is based on random 
assignment [24] which can be described using a probabilistic mapping device which randomly assigns bits of the 
transmitted codeword to M blocks. More precisely, the random mapper takes a bit of the codeword C and sends 
it to slot j with a predetermined assignment rate r, . Bit assignments are independent and, thus, the expected (and 
asymptotic as ^ oo) number of bits assigned to block j is TjN. 

After random mapping, each block of bits forms a codeword obtained by puncturing. Let Cj denote the punctured 
codeword corresponding to slot j. For a given good (mother) code ensemble [C] and the (random) assignment rate 
Tj, Theorem [2] shows that the punctured code ensemble exhibits the SNR threshold behavior for an AWGN channel, 
i.e., if the received SNR is larger than the punctured code threshold x('^j)' the average ML decoding word error 
probability decays to zero as the codeword length N tends to infinity. As illustrated in (|2T]) . the punctured code 
threshold xi^j) can easily be calculated based on the mother code threshold c^. Hence, this threshold behavior 
allows for adaptively scheduling the cooperation in two stages described below. 




o 



o 



o 

□ 

destination 

O 



collaborating cluster 

a. slot 0: sender broadcasts message 
Fig. 5. Cooperative coding scheme 



sender 




O 

H 

destination 



O 



collaborating cluster 

b. slot 1 - slot M — 1: reliable nodes relay information 



A. Broadcast Stage 

As shown in Fig. |5]a, the sender broadcasts its information by transmitting Co during slot 0. Each helper listens 
and attempts to decode this message. More precisely, helpers estimate the instantaneous received SNR from the 
incoming signal and compare the SNR with the punctured code threshold x(''"o) corresponding to Cq at the beginning 
of the broadcast stage. Let C 5 denote the reliable set of cluster members whose sender-to-helper instantaneous 
SNR Oqj > x{'^o)- The element of F is referred to as the reliable node. The punctured code threshold theorem and 
quasi-static channel assumption imply that reliable nodes can be guaranteed to decode the message successfully 
under ML decoding. Next, each reliable node sends an ACK back to the sender over a fast and error-free feedback 
channel without a need for prior decoding. During this slot, reliable nodes listen and decode the received signal. 
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TABLE I 

Dynamic schedule {M = 5 and ^ {1,3}) 



Slot 





1 


2 


3 


4 


Transmitted block 


Co 


Ci 






C4 


Transmitter 


Node 


Node 1 


Node 


Node 3 


Node 



B. Forwarding Stage 

Fig. |5]b illustrates the transmission interval corresponding to slots 1 through M — 1 (here, dark solid circles 
represent reliable nodes). Node k £ T, re-encodes and partitions received information in the same manner as done 
by the sender and, consequently, relays Ck (corresponding to block k) to the destination in slot k. The sender 
transmits the left-over blocks in the remaining transmission slots based on received ACKs. The signal received at 
the destination corresponds to an IR scheme with a fixed number of retransmissions, where a retransmission may 
experience a different channel quality. 

Example 3 Let's assume M = 5, Node is the sender, the collaborative cluster S = {1, 2, 3, 4}, and the reliable 
set T = {1,3}. Based on the ACK-based feedback from the cluster, the dynamic schedule is shown in Table H 
Since Nodes 1 and 3 are reliable nodes, they send Ci and C3 at Slots 1 and 3, respectively. The remaining blocks 
are transmitted by Node at Slots 0, 2, and 4. 

Remark 3 Decoding Delay: We note that each reliable node needs to decode the message before it relays this 
information. This results in a decoding delay in the transmission. If one of the reliable nodes is scheduled to send 
the message in slot 1, it requires some extra time td ■ T between slot and slot 1 due to a decoding latency. In 
Appendix IB we propose both early stopping and a threshold adjusting technique to reduce this time in practice. 
For simplicity of the error performance analysis, we assume r/5 = in Sec. |Vl 

Remark 4 Decoding Failure: Since the assignment rate tq is predetermined, each helper can pre-calculate the 
punctured code threshold xi^o) by using ((2T]) . On the other hand. Node j may estimate the instantaneous received 
SNR Oqj at the beginning of slot 0. In our proposed scheme, if Oqj > x{to)^ Node j is called reliable node and 
sends an ACK message to the sender. Next, if Node j fails in decoding the message, this node will stay silent 
during the forwarding stage. We refer to this event as the decoding failure. However, Theorem [2] implies that the 
probability of such events approaches zero as the codeword length N 00. In this paper, we focus on long codes. 
Hence, we neglect the decoding failure event in the error performance analysis. 
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V. IR Cooperative Coding Performance Based on Code Outage 

The decoding is performed at the destination upon completion of M transmission slots. The proposed cooperative 
coding scheme implies that the received signal is always the mother codeword C modified by the fading 
channel. Moreover, all communication links experience independent quasi-static Rayleigh fading channels. Thus, 
the codeword C is, equivalently, transmitted over M slots and experiences |jr| + 1 independent channel gains. 
Consequently, we study the performance of codes transmitted over a block fading channel [15], [20]. 

A. Code Outage for a Q-Block Fading Channel 

Here, we consider the block-fading channel model [20] with Q fading blocks (a group of Q blocks will be 
referred to as a frame), where the fading coefficient is essentially invariant during a single block and different 
from one block to another. Let SNRj and Uj be the the average received SNR and the channel power of block j, 
respectively. The Bhattacharyya parameter 7^ is a function of SNRj and Uj, i.e., 

7j = exp(— i^j • SNRj) for j = 0, . . . , Q — 1. 

Hence, the average Bhattacharyya parameter over Q blocks 

Q-l Q-1 

is a function of the random vector v = {uq, . . . , t'g-i} and, thus, for a given good code, there is a non-negligible 
probability that the effective Bhattacharyya distance — In 7(1/) is less than the code threshold c^, termed code 
outage probability. Thus, the average frame error probability is a function of both the fading distribution and the 
threshold of the code ensemble [C]. More precisely, the average ML decoding word error probabiUty for a good 
code ensemble [C] transmitted over a Q-block fading channel can be bounded as follows: 

= P{error(A^), — In 7(1^) < c^} + P{error(A^), — ln7(i>') > c*} 

< P{-ln7(iv) < c^} + P{error(iV) | - In 7(1/) > c^,} (24) 

where error (A^) stands for the event of the decoding error for the code of length A^. Theorem [1] implies that the 
second term of (l24l ) approaches zero as the code length N increases. Hence, the following error rate bound holds 

Pwil, N) < F{^{u) > exp(-c.)} + 0(1) (25) 
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/ N TV 

where o(l) — > 0. We focus on long codes in the following analysis. To simplify notation, we will omit N and 

0(1). 

B. IR Cooperative Coding Performance 

Here we study the FER performance of the IR cooperative coding scheme for a quasi-static frequency-fiat 
Rayleigh-fading channel based on the code outage upper bound ( |25] ). 

In slot 0, the sender (Node 0) broadcasts its information by sending the punctured codeword Cq. The channel 
powers vqj, j = 1, . . . , M — 1, are i.i.d. exponential random variables invariant during each transmission period. 
The reliable set is now randomly distributed over the collection of 2^^~^ subsets of S with probability 

P(.^) = n > ^(^o)} n Pi^o,, < x(ro)} 

= n«^P[-^(^o)SNRoj] n {l-exp[-x(To)SNR-]]}, (26) 
jar j&r- 

where = S\T. 

For a given J^, the IR cooperation scheme allows the mother codeword to be transmitted to the destination 
(Node M) over M slots with |JP"| + 1 independent quasi-static fading gains. Hence, the slot i Bhattacharyya 
parameter is 

exp(-6'i,M) i G 

li = { (27) 
exp(-eo,Af) iG^^U{0}. 

Consequently, the Bhattacharyya parameter averaged over M slots is now 

7(1^, JF) = (1 - ^ Ti) exp(-6'o,Af) + ^ Tiexp(-6'i,M) 

= (1 - ^ Tj) exp(-z^o,A/SNRo,Af) + exp(-fj,MSNRi,M), (28) 

where v = {uq^m, ^i,m ■ ■ ■ , i^m~i,m} is a random vector M-tuple with an independent exponential distribution. 
The bound (1251 ) implies that the conditional average word error probability given a reliable set can be bounded 
as 

» M-l 

Pw{l I ^) < P{l(.J^,^) > exp(-c.)} = / n e-"''^' du 

■'^ i=0 

= g{M, J^,SMR) (29) 

where A ={v : 7(1/, JP") > exp(-c^)}, SNR = {SNRo,Af , . . . , SNRm-i,m}, and C/(M,J^, SNR) is referred to 
as the code outage probability for a given reliable set J^. The ML decoding FER for cooperative coding scheme 
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averaged over all possible reliable sets can be bounded as follows 

FerW= ^ p(^)p^(^|^)< P(^) .g(M,^,SNR) 

all possible all possible 

where the superscript (M) represents the number of (potential) transmitting nodes. 



(30) 



Example 4 (M = 1) The case M = 1 is equivalent to the traditional direct transmission between the sender and 
the destination, and the single-hop FER is 

FERW <g(l,0,SNRo,i) 
= P{z.o,iSNRo,i < 

= l-exp(-c*SNRo}). (31) 

Example 5 (M = 2) Here we consider the IR cooperative coding scheme with a single helper. In this case, the 
scheme FER is 



FER(2) <F{T = $)Pw (7 I 0) + P(^ = {l})Pw (7 I {1}) 



1 - e 



where (assuming ro,ri < exp(— c^)) 



x(r)SNR-; 



g{2, 0, SNR) + e- xMSNRoj g^2, {1}, SNR), 



g{2, 0, SNR) = 1 - exp(-c^SNRo j 
g{2,{l},SNR) = l-ij- 



^ -exp(-c^) - ro2;SNRo,2 ^ i/snRi,. 



exp[— x(to) • SNRq2], and the intermediate steps for deriving (l34l) are given in Appendix] 



(32) 



(33) 
(34) 



In general, (1301 ) cannot be calculated in a closed form and one needs to resort to numerical integration methods. 



VI. Asymptotic Analysis 

In this section we consider several different cooperation scenarios and derive asymptotic (in SNR) FER bounds, 
which have a closed form. For simplicity, we assume 



To,..., TM-i > 1 - exp(-c^). 



(35) 



i.e., each randomly punctured code is self-decodable (with probability one). Next, we refer to r = 
max[do,i, • • • ,do,A/-i] as the sender-to-cluster distance, D = do^M as the sender-to-destination distance, and 
d = max[di_Af , • • • ,dM-i,M] as the cluster-to-destination distance as shown in Fig [T] Similarly, we define the 
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sender-to-cluster SNR, the sender-to-destination SNR, and the cluster-to-destination SNR as 

p^E-r-^, r] = E-D-^, and X = E-d-^. (36) 

Computing FER bound ( [30l ) requires integration in ^(M, JT, SNR), which presents the code outage probability 
for a given J^. The following theorem which is the basis of our asymptotic analysis, among their contributions, 
helps in avoiding this integration. 



Theorem 3 Consider Q independent random variables (pi,. . . ,(pQ with the following properties: 



< (/>m < 1 for m = 1, . . . , Q 



where the probability distribution of 4>rn is a function of Am such that 



lim Xm P [4>m > c] = - In c 



and < c < 1. If Ti, tq > 1 — c and J2m=i '^m = 1, then 



Q 



lim sup AmP 

Ai,...,Aq^oo^^-^ 



E 



Tm.0m. ^ C 



1 T 



m=l 



(1 - Tm) ' 



where Ai, . . . , Aq ^ oo means Ai oo, . . . , Ag ^ oo. 

Proof: The proof is provided in Appendix |H] based on the induction method. 



(37) 



(38) 



A. Transmitter Clustering 

In the transmitter clustering scenario we assume that M — 1 helpers are very close to the sender so that r ^ 0. 
For this setting, we call the sender-to-cluster channel ^/fy reliable in the sense of Y'{J^ = {1, . . . , M — 1}) = 1, 
i.e., all helpers are reliable nodes with probability one. Thus, the cooperation scheme FER can be written as 

FER?"'^ =P^(7|.F = {1,...,M-1}) 

<P{-ln7(iy,{l,...,M-l}) <c.} 

= g{M, = {1, . . . , M - 1}, SNR = {EdTj;,]). (39) 
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Now, let's consider the large SNR case. Note that the exponential distribution of f j m implies that 



lim SNRiA/P 

SNRi,Af^oo 



lim 



exp(-i/j,MSNRi,M) > exp(-c^ 
l-exp(-c.SNR-]^) 



l,M 

=c^. (40) 
Thus, (ED, (Hill, (l35]l, and Theorem |3] imply 

M-l M-1 

limsup n SNR,,Af •a(Mai,...,M-l},SNR) < — n x(Ti). (41) 
SMR-c» tl, to 



For large enough E, we can rewrite (1391 ) as 

M-l 

FER?^) XfT n ^(^OSNRrl, (42) 



M! 

i=0 



M-l 



■ i=0 



(43) 



= ^M]^ n (44) 

■ i=0 

where means that the inequality holds for sufficiently large E, and the last step is based on the triangle inequality 

D — r < di^M < D + r and r ^ 0. 

B. Receiver Clustering 

In the receiver clustering scenario we assume that M — 1 cluster members are very close to the destination so 
that d 0. Note that (l35l) implies that each block (punctured codeword) is self-decodable with probability one. 
Hence, the code outage probability is zero for any nonempty rehable set J^, i.e, 

g(M,J^,SNR) = for J^/ 0. (45) 

Therefore, we can bound the cooperation scheme FER as 

FER^''') = P{T = 0)p[^l (7 1-^ = 0) 

< P(.F = 0)g(M,.F = 0,SNR) 

M-l 

= n {1 - exp[-x(To)SNRo j] } [l - exp(-c.r?-i)] 
i=i 

M-l 

= []{!- exp[-x(To)i?-^4J } [1 - expi-c.E-^D'^)] (46) 



Again, we focus on the large SNR case. Note that 



lim SNR- ri-exp(-aSNR-^)l =a for a > 0. 

SNR^oo '- 

Hence, for large enough E, we can rewrite (1461 ) as 

M 

FERf) <E [x(ro)]*'-'c..nSNRo,] 

_ [x(ro)]- 



\M-1 M 

C- 



[Xiro)] 



EM 

M~l 



n 4. 



2^ _ jjML 



EM 

where the last step follows from the geometric property D — d < doj < D + d and d 0. 



(4 



(4 



C. Cluster Hopping 




sender destination 
collaborating cluster 

Fig. 6. Cluster hopping 

Here, we assume that (io,i = • • • = do,A/-i = r > and di^M = • • • = dM-i,M = d > 0. The distances d, 
and D satisfy the triangle inequality as shown in Fig. |6] By using (l26l l in the high SNR regime, we have 



lim TT SNRo, pm = lim /^"(i-^i+i)pm 



lim I exp[-x(ro)/5 

[^(^„)]M-(|^|+1). 



-11 I J 1 - exp(-x(To)/3 ) I 



(4 



Moreover, dHJ, dUl, (l35]l, (gOll, and Theorem |3] imply 



limsup JJ SNR,,M-a(M,J^,SNR) =limsup(Al-^lr/)-g(M,J^,SNR) 

SNR-.00 ^^^^^Q^ ' ^ " 



< 



A, T]-^00 



m + 



(5 
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Note that 



SNRo^M = V = ED-^, SNRoj = p = Er-^, SNR^- m = A = -Efi"^, j = 1,...,M -1. 



Hence, for large enough E, we can rewrite (1301 ) as 

FER<") Y.[ ''''"°""'"ij!;^|"+'i)r n x(r.)|(pl^l+'-"A-l^l ,-1) (51) 



geometric distance profile 



coding advantage 



D. Diversity Gain 

Following [27], the diversity gain is defined as 

divA ^1^^. (53) 

SNR^oo logSNR 

Since our collaborating model is a distributed multiple-input single-output (MISO) system, the maximum achievable 
diversity gain is M. Equations (l42l ). WT\ . and (ISTl ) illustrate that all of the three discussed scenarios: transmitter 
clustering, receiver clustering, and cluster hopping can achieve the full diversity gain, i.e., div = M, in high SNR 
regime. 



E. Cooperative Coding Gain 

For small c^,, we can build the following simple relationship between the punctured code threshold and the mother 
code threshold c^,. Equation (|2T] ) implies 

lim ^ = lim (c*)-i In ^— = -. (54) 

c*^o c»^o exp(— c*) — (1 — r) r 

Thus, we can rewrite (l52l ) as 

FERI") E{ M-,|:r|..,., ^ — — } ' (r^'-'^'-'d'^' ■ E'-. (55) 

(1 - Eie.F^i)(nie.F7-i)(|.F| + 1)! J 

where <e,c^, means that the inequality holds for sufficiently large E and small c*. 
Example 6 (M = 1 limiting case) 

FER(i) <ij - = ^D^. (56) 

?7 E 
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Example 7 (M = 2 limiting case) 



TqTJP 2tqTi 7]X \E 

Similarly, (|44l) and (1481 ) can be rewritten as 



To 2toTi 



(57) 



M 



-(M-1) I c^D 



fER'^" <E,c. r, [::^] . (59) 



In [28], the author defines the cooperative coding gain as 



COp= lim -..p ■ (60) 

SNR^oo SNR 

Let the sender-to-destination SNR be the basis, i.e, r/ = SNR. Bounds (1581 ). ( [59l ). and (1551) imply that the coding 
gains of transmitter clustering, receiver clustering, and cluster hopping schemes satisfy 

Urn • COp^^^) > M! J] rj , (61) 
Umc..COpf)>rr^-^)/^, (62) 

Inequality (|63] ) illustrates that, in general, the cooperative coding gain is a function of the code parameter c^, the 
cooperation scheme parameters {tj}, and the geometric distance profile (r, d, D) of the network. 

VII. Simulations and Discussions 

A. IR Cooperative Turbo Coding 

In this subsection we study the error performance of the IR cooperation scheme based on the turbo code described 
in Example [T] PER simulations consider binary antipodal signaling and an independent flat quasi-static Rayleigh 
fading for each link. Each receiver has perfect channel state information and employs coherent detection. All 
receivers employ a multiple turbo decoder based on the triangle iterative decoding algorithm [29]. 

Here we consider a M = 5 collaborative network and assume SNRo,i = • • • = SNRo,4 = p and rj = SNRo,5 = 
• • • = SNR4 5 = A. Thus, the FER performance of cooperative turbo codes is a function of both cluster-to-destination 
SNR p and sender-to-cluster SNR A. Fig. |7] depicts the FER for p = —15, 0, 15 dB and as a function of A from 
—2 to 16 dB. On the other hand, in Fig. [H we fix A = 2, 6 dB and study the FER performance vs. sender-to- 
cluster SNR p. For these two cases, we compare the simulation result with the analytic upper bound (l30l ) and the 



p£p-l/diV 
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o simulation (cooperative) 

— upper bound 

— asymptotic bound 
* simulation (direct) 



cc 

LU 




-2 2 4 6 8 10 12 14 16 
cluster-to-destination SNR A, (dB) 

Fig. 7. FER vs. A (M — b, p — —15, 0, 15, and mother turbo code of rate R — 1/7 and length N — 5376) 
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sender-to-cluster SNR p (dB) 



15 



Fig. 8. FER performance vs. p (M = 5, A = 2, 6, and mother turbo code of rate i? = 1/7 and length N = 5376) 



asymptotic bound (|55] ). Figs. [V] and [8] also depict the simulation result of the direct transmission as a benchmark. 
We observe that the upper bound (l30l ) accurately predicts the cooperative coding performance and the asymptotic 
bound (22]) converges to the bound (l30l ) for medium and high SNR. This observation enables us to estimate the 
FER performance as a function of p and A by combining (l30l ) and (1551 ) in Fig. |9l where we use the bound (l30l ) for 
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low SNR and employ the bound ( [55] ) to simplify the computation for medium and high SNR. 
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sender-to-cluster SNR (dB) 



5 



Fig. 9. FER performance as a function of p and A {M ~ 5 and mother turbo code of rate R~ 1/7) 

B. Collaborative Cluster Size 

Here, we study the effect of the collaborative cluster size on the FER performance of the transmitter clustering 
(the sender-to-cluster distance r ^ in this scenario). We assume that nodes have limited battery energy. In this 
case, achieving high transmission energy efficiency is more important than maximizing the diversity gain. Our 
approach is to assume that the allowable FER is e, which guarantees the quality of service (QoS), and to determine 
the e-achievable transmission energy by applying the asymptotic bound studied in Sec. |Vl] The closed form bound 
predicts well the error performance for medium and high SNR. 

Let To = • • • = TM-i = now, (1581 ) implies 



'strictly speaking, tlie bounds l |55t . l |58l l, and l |59l l are based on the large SNR assumption. However, through simulations, we observe 
that the asymptotic bounds also works well for the medium SNR. On the other hand, these asymptotic bounds can be expressed in a closed 
form, whereas, the calculation of the bound l |30l l requires numerical integration method. Thus, here and hereafter, we use these asymptotic 
bound to estimate the FER performance. 




(64) 



where E'(*^) is the e-achievable transmission energy. To satisfy the QoS requirement, we requirq^ 




(65) 
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TABLE II 
Energy saving vs. M 



M 


2 3 4 5 


FIRF 


Estimated C/^*^) (dB) 


8.4 11.1 12.4 13.2 


< 20 



Based on Stirling's approximation, the e-achievable transmission energy is 

E(M) ^ 

To illustrate how much energy can be saved using the IR cooperative transmission, we consider both direct 
transmission and transmission over a fully interleaved Rayleigh fading channel cases as benchmarks. For direct 
transmission (i.e, M = 1), by using the bound (l56l) in Example [6l the e-achievable energy is given by 

^(1) ^ (67) 
e 

For a fully interleaved Rayleigh fading channel, the Bhattacharyya parameter is a function of the sender-to- 
destination SNR rj (see [17] for the detail) as follows 

7(™^) = (68) 

where (FIRF) stands for fully interleaved Rayleigh fading. Theorem [T] implies that the asymptotic word error rate 
approaches zero (as N oo) if 7(^^^^) < exp(— c^,), i.e., the sender-to-destination SNR 

■q > exp(c^) - 1. 

Note that rj = E ■ D^^. Thus, for this channel, the reliable transmission energy threshold is defined by 

^(FiRF) A [exp(cj - 1] • D^. (69) 

Let = denote the transmission energy saving. Equations (l66l ) and (l67l ) lead to 

which illustrates the fact that the energy saving C/(^) is a function of only M and e for the transmitter clustering 
scenario and does not depend on the code threshold Cj, and the sender-to-destination distance D. In other words, the 
energy saving of the IR cooperative coding scheme is universal for all good code families and sender-to-destination 
distances. Next, we set e = 0.01 and numerically compute the estimated energy saving in Table |ll] based on the 
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Fig. 10. FER vs. transmission energy in the transmitter clustering scenario (D=l, P{\J- \ = M — 1) = 1, and mother turbo 
code of rate and length N = 5376) 



approximation (TTOl) . We compare the calculation result with the fuUy interleaved fading channel savings 

jjiFIRF) ^ 

^(FIRF) [exp(c.) - 1] • e ^ / ' ^ ^ 

where the inequality follows from exTp{x) > 1 + j; for x > 0. Fig. [10] illustrates the simulated FER performance 
versus transmission energy E as well as the upper bound (l30l) and its asymptotic version (1551) for D = 1 and 
r ^ 0. In Fig. \T0\ we also compare the FER performance of cooperative transmissions vs. transmission over a fully 
interleaved Rayleigh fading channel. In the latter case, the error performance exhibits a threshold behavior and the 
reliable transmission energy is described in ( [69l ). We observe that the energy saving obtained through simulation (in 
Fig. [TOb and the estimated J7(*^) (in Table HIl) illustrate an excellent match. Furthermore, both Table HIl and Fig. [TOl 
illustrate the fact that, although the cumulative energy saving increases with M, the rate of increase drops quickly. 



C. Normalized Cluster to Destination Distance 

Here, we assume cLq i, . . . , do^M-i = di^M = • • • = dM-i,M = d, D k, r + d and ri = • • • = tm~i- We move 
the collaborative cluster from the sender towards the destination, and evaluate the energy saving in terms of the 
normaUzed cluster to destination distance k = r/D. Here, we use the similar approach as one used in the pervious 
subsection. We assume that the allowable FER is e and calculate the the e-achievable energy E^^'^^ based on the 
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asymptotic FER bound dSSl ). i.e., 



Hence, we have 



where rj- = (1 — to)/(M — 1). Next, the transmission energy saving is given by 



-l/M 



= = ) ,M-i V ^- ^ ' \ (74) 



_ J ^^iZ,^ /Af - 1\ ^L(Af-fc -l)(i _ ^)Lfc 

2^ I A; j Af-(fc+i) 

Equation (l74l ) illustrates that does not depend on the code threshold and the sender-to-destination distance 

D for a given k. In this sense, we claim that the transmission energy saving of IR cooperative coding schemes is 
universal. 

Based on (1741 ). we depict the transmission energy saving U^'^^^ in Fig. [TT]as a function of the normalized cluster 
distance k for M = 2, 3, 4, 5, 10, where e = 0.01, path loss exponent L = 3, and < k < 1. Here, we consider 
two cases: fixing tq = l/M and choosing optimum tq for a given k. In Fig. \TT\a, we observe that U^^^^ is below 
f/^^) when K — > 1. The reason is that here the assignment rate is fixed to tq = l/M. We note that k — > 1 implies 
that the collaborative cluster is close to the destination and far from the sender. In this case, the transmission in 
slot (i.e., broadcast stage) is more important (see the receiver clustering case in Sec. IVI-BI and the asymptotic 
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bound (|59l)). and larger tq helps in increasing the expected number of reliable nodes which relay the message to 
the destination at the forwarding stage. Hence, we obtain a better energy saving performance for tq = 1/4 than for 
To = 1/10 when /t ^ 1. As a comparison, we depict the energy saving by optimizing tq for a given k in Fig. [TTlb. 

VIII. Conclusion 

In this paper, we study a cooperative coding scheme for multi-hop wireless networks with incremental redundancy. 
Threshold analysis of good code ensembles enables the analysis of the IR cooperation coding scheme. 

First, we study the threshold behavior of good code ensembles for BISOM Channels. A general relationship 
among the channel capacity, the Bhattacharyya rate, and cutoff rate is established for BISOM channels. Based on 
this relationship and the modified Shulman-Feder reliable channel region [15], a simple code threshold is proposed. 
The reliable communication constraint based on the simple code threshold has the same form as the UB code 
threshold introduced in [17], but the former result is tighter by almost 1 dB for an example of a rate 1/7 turbo 
code on an AWGN channel. 

Next, we consider the fading channel. Based on the code outage concept, a simple threshold bound on FER and 
its asymptotic version are presented for the IR cooperative coding scheme. The asymptotic bound is a closed form 
function of the simple code threshold, average channel SNR, geometric distances, and the size of the collaborative 
network. It enables analysis of the diversity and coding gains for three different scenarios: transmitter clustering, 
receiver clustering, and cluster hopping. 

Finally, we simulate the IR cooperative turbo code performance based on the iterative decoding. Remarkably, 
our analytical FER bound accurately predicts the simulated cooperative coding performance at any SNR and the 
asymptotic bound agrees very well with the simulation results at medium and high SNRs. We further discuss 
the transmission energy gain in terms of collaborative network size and normalized distance for the transmitter 
clustering and cluster hopping scenarios, respectively. In both cases, we observed that the energy saving does not 
depend on the sender-to-destination distance and the code threshold. In this sense, we claim that the energy saving 
of IR cooperative coding scheme is universal for all good code families and all initial non-cooperative hop-distance 
selections. 

Appendix 

A. Modified Shulman-Feder Reliable Channel Region 

In this appendix, we state for reference the Modified Shulman-Feder (MSF) reliable channel region for good 
binary code ensembles transmitted over a set of Q (independent) parallel channels (see [15] for the detail). 
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Theorem 4 (cf. [15, Theorem 5]) Consider a good binary linear code ensemble [C] of rate R whose transmission 
takes over Q BISOM parallel channels. Assume that the coded symbols are randomly and independently assigned 
to these channels, such that each bit is transmitted over Channel j with a-priori probability tj, for j = 0, . . . , Q — 1, 
where tj > referred to assignment rate and Y^'j=Q tj = 1- Let Cj and jj be the channel capacity and the 
Bhattacharyya parameter of Channel j, and cp and be the restriction UB code threshold and the SF distance 
defined in (fT6l ) and (ITtI ). respectively. Then, if the average channel capacity and Bhattacharyya parameter over the 
Q parallel channels satisfy 

Q-l Q-1 

Tj'jj < exp(— Cp) and ^ TjCj > R + ^p, (75) 
j=0 j=0 

the average ML decoding word error probability decays to zero as the codeword length approaches infinity. 

Different from the condition ([T9l ) for the simple code threshold, the above MSF reliable channel region (175] ) 
requires satisfying two constraints associated with two pairs of channel/code parameters. Hence, this calculation is 
complex, in particular, when the condition (175] ) is employed with some practical communication schemes. 

B. Some Useful Inequalities 

Proposition 1 l + 6^<2''<l + 6forO<6<l, where the two equalities hold simultaneously when 6 = 
and 6=1. 

Proof: Let g{b) = 1 + 6^ — 2''. We check the derivatives of g{b) in the region 6 G [0, 1] as follows 
^'(6) = 26-ln2•2^ 

g"{b) = 2 - (ln2)2 . 2'' > 2 • [1 - (ln2)2] > 0, for < 6 < L 

This means that g{b) is strictly convex. Note that g{0) = g{l) = 0. Thus, g{b) < for 6 G (0, 1) and the first 
inequality of Proposition [T] holds. 
Next, let /(6) = 1 + 6 - 2\ we have 

/'(6) = l-ln2•2^ 

/"(6) = -(ln2)2 -2* < forO<6<L 

This implies that /(6) is strictly concave. Since /(O) = /(I) = 0, we have /(6) > for 6 G (0, 1) and the second 
inequality of Proposition [T] holds. ■ 

Proposition 2 2~" + 2~^/°- < 1 for a > 0, where the equality holds when a = 0, 1, and oo. 
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Proof: Consider /(a) = 2~" + 2^^/". Due to the symmetry property of the function /(•), i.e., /(a) = /(1/a), 
we only need to prove /(a) < 1 for < a < 1. 

Note that /(a) is differentiable and bounded (0 < /(a) < 2). Thus, the maximum value of /(a) is either at a 
boundary point or at a stationary point. Now we check all such kind of points. For boundary points, clearly, 

/(a) = 1 when a = 0, 1. (76) 

By the definition, the stationary point s satisfies 

/'(a) 1^^^ = - ln2 • (2-^ - ^2-^/^) = 0. (77) 

This implies 2^^/* = ■ 2^^ . Thus, for a arbitrary stationary point s, we have 

/(s) = 2-'' + 2'^/^ = 2~^ + 5^ • 2-^^ = 2~*(1 + s^) < 1 forO<s<l (78) 

where the inequality follows from Proposition [T] Hence, (1761) and (ITSl) imply the desired result. ■ 



C. Proof of LemmaU} 

The proof of Lemma [U is based on Jensen's inequality [30] and Proposition [U |2] in Appendix |B] 
Proof: First, we prove B < C. Note that 

C-i? = -E|log2 

Following Jensen's inequality [30], we have 



p{Y\0) +p{-Y\0) 2-Vp(-^|0)/p(y|0) 
p{Y\0) 



C-B> -log2 |^[p(y|0) +p(-2/|0)] • 2-Vp(-2/|0)M2^|0)1 



logs I Piy = 0|0) + ^ [p{y\0) + p{-y\0)] 

[ y>0 



Vp(-j/|o)/p(s/|o) _^ 2"Vp(j/|o)/p(-J/|o) 



By using Proposition |2l we have 



(79) 



y/pi-y\0)/p{y\0) _^ 2-VpW7?K^^ 



< 1 
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Hence, dTOb can be bounded as 



C-B> -\og^[p{y = 0|0) + j;[p(?/|0) +p(-y|0)]} 



y>0 
= - l0g2 1 

= 0. (80) 
Now, we prove B > Rq. Following Proposition [H we have 

7<log2(7 + l) forO<7<l, 
By the definition of the cutoff rate in we can bound the Bhattacharyya rate as 

S = l-7>l-log2(7 + l) = ^o. 
By combing (|80l ). we have the desired result. ■ 



D. Proof of Lemma^ 

The proof of Lemma |2] is by contradiction. 

Proof: We assume that (fT3l) does not hold. Then, there exists a positive eo such that 



> 1 - exp(-co) + eo. (81) 

Let's consider a binary erasure channel (BEC) with erasure probability p = exp(— cq) — eo. Then, the channel 
capacity and the Bhattacharyya parameter are 

C{p) = I - p = I - exp(-co) + eo and j{p) =p = exp(-co) - eo- 

Since eo > 0, the UB reliable communication condition (fT4l ) is satisfied, i.e., ^{p) < exp(— co). Hence, the decoding 
error probability approaches as ^ 0. Now, the converse to Shannon's channel coding theorem [30] implies 

R < C{p) = 1 - exp(-co) + eo, 

which contradicts (ISTTl. ■ 



E. Proof of Theorem [7] 

For the proof of Theorem [T] we proceed in the following two steps. We first show the existence of by using 
Lemma |2l Next, we prove the main result based on the MSF reliable channel region theorem (see Appendix lAl ) 
and Lemma [T] 
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Proof: Let 

O = {P : 1 - exp(-cp) > i? + and < P < 0.5}. 

To prove the existence of c^, we need to show that the set is not empty. Let's consider a particular weight 
partition that P = 0. In this case, Equations (fT6l ) and (fTTl ) imply cp = cq and = 0. Moreover, Lemma |2] leads 
to 

1 - exp(-cp) >R + ip for P = 0. 

Hence, the set 7^ and is well defined. 

Next, we prove that, for V P G Vt, Pw{l, iV) ^ if 

7 < exp(— Cp) . (82) 

The condition (l82l ) implies 

1 -7 > 1 -exp(-cp) > P + Cp forPGfi 

where the second inquilinity follows from the definition of Q.. On the other hand. Lemma [T] implies that the average 
channel capacity 

Q-l Q-1 
C^5;T,C7,>^(l-r,7,) = l-7. 

j=0 j=l 

Therefore, 

C > P + ^P. 

By combining (l82l) and Theorem IH we have the desired result Piy(7, A^) 0. ■ 

F. Discussion on Reducing Decoding Delay 

Here, we propose "early stopping" and "threshold adjusting", which in practice neutralize the effect of decoding 
delay. 

Without of loss generality, we assume that Node 1 is a reliable node and scheduled to send the message in 
slot 1. The key of early stopping is that Node 1 does not need to listen to the message during the whole slot 
[0, tqT], instead. Node 1 can stop listening early and begin to decode after it receives enough information. Let 6*0,1 
be the sender-to-helper instantaneous SNR of Node 1. By the reliable node definition (see Sec. IIV-AI ). we have 
^0,1 > x(''"o)- Hence there exists a effective listening period r' satisfying 

^0,1 > X(r ) > x(ro). (83) 
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Equations ( [83] ) and (|2TI ) imply that r' satisfies 



1 - exp{-c^) , 

<T < Tq. 



1 - exp(-6'o,i) 

Node 1 begins to decode after receiving t'N bits. Theorem |2] implies that the "early stopping" rule can guarantee 
ML decoding successful at Node 1. Now, the required extra time is td + t' — tq. Clearly, the proposed strategy 
helps in reducing decoding delay latency. 

When TD <^ To, an alternative method is threshold adjusting. Here, we reset the punctured code threshold to be 
x{to—td)- We say that Node 1 is a reliable node (scheduled to send the message in slot 1) when the sender-to-helper 
instantaneous SNR 

Oo,i>x{To-rD). (84) 

Node 1 listens to the message during the slot [0, (tq — T£))T] and begins to decode after receiving (tq — T£i)N 
bits. Now, there is still t^T in slot for decoding process. 



G. Proof of Equation f l34D 

Here, we provide the intermediate steps for deriving (l34l) . 
Proof: Let 

X = exp(-fo,2) and y = exp(-z^i^2)- 

Recall that fading power gains i/i 2 and 1^0,2 are independent exponentially distributed with mean 1. Hence random 
variables x and y are independent uniformly distributed over (0, 1]. By the definition ( |29l ). we have 

g{2, {1}, SNR) = P{ro7o + ri7i > exp(-cO} 

= 1 - P{ro7o + ri7i < exp(-c^)} 

= 1 - P{ro7o + Ti < exp(-c*)} - P{ro7o + n > exp(-c^), ro7o + ri7i < exp(-c^)} 

= 1 - P{7o < exp[-x(ro)]} - p{7o > exp[-x(ro)], 71 < ^^Eizf^l^^} (85) 

where the third equality is due to < 71 < 1, and the last equality is based on the definition of punctured code 
threshold (EB- Note that 

70 = exp(-i/o,2SNRo,2) = x^n^''^^ 

71 =exp(-zyi,2SNRi,2) = y^^'^^^ 
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Hence, the second term of ( [85]) is given by 



P{7o < exp[-x(To)]} = P{x < u} 



to 



(86) 



where oj = exp[— x(to) • SNRq 2] • Next, the third term of (1851 ) can be rewritten as 



p|7o > exp[-x(ro)], 71 < 



exp(-c^) - To7o 



F < X > u), y < 



exp(-c*) - tqx 



SNRo- 



Tl 



1/SNRi 



exp(-c^) - tqx 

Tl 



SNRo,: 



1/SNRi 



dx 



Combining (1851 ). (f86b . and (f87b together, we have the desired result ([34. 



(87) 



H. Proof of Theorem \3\ 

The intuition of the proof of Theorem [3] is described as follows. Since Q is a positive integer, we derive the 
result based on an induction method. First, we consider the Q = 2 case and derive Lemma [3] based on a more 
general setting. Next, we prove the main result of Theorem [3] We show that if the hypothesis (1381) holds for the 
integer Q = j — 1, it is true for the next greater value Q = j. 

Lemma 3 Consider two independent random variables and (1)2, where < (pm < 1 and the probability 
distribution of is a function of A^, for m = 1, 2. Assume that 



lim AiP[(/)i > c] = /(c), limsup A2P[(?!)2 > c] < h{c), 
where /(c) and h{c) are monotone decreasing and integrable, and /'(c) is integrable. If 



(88) 



< r, 1 - r < c, 



then 



limsup AiA2P[r(/)i + (l-T)02 >c] < - / h(^:^^) f'iz)dz 

Ai,A2^oo Jq V 1 - r 



(89) 



where q = [c — {I — t)]/t. 

Proof: [Lemma 121 The proof follows the approach in [3]. First, let Z = {zo,...,zl} for some finite L 
be any partition of the interval [q, 1] with zq = q and zl = I. Next we obtain an outer bound on the event 

{T(j)i + (1 - r)(/)2 > c} as 



i=l{ ^ . 



(90) 
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Since (pi and (/)2 are independent, we have 



Zi^l <(t>l< Zi, 4>2 > 



C ~ TZj 
l-T 



P[01 > Zi-l] -P[</>1 > Z,]}P 



h > 



C - TZj 

1 -r 



(91) 



Hence, 



limsup AiA2P[t</)i + (1 — T)(p2 > c] < limsup A1A2P 



i=l 



Ai ,A2— ►oo 



Zi-1 < (pl < Zi, (j)2 > 



C - TZj 
l-T 



Li 

< limsup Xi{f[4>i > Zi^i] - P[0i > zA ■ A2P 



t>2 > 



C - TZj 
l-T 



C - TZi 



(92) 



Note that ( |92l ) holds for all partitions Z of the interval [q, 1]; and /(c), /'(c), and /i(c) are all integrable, the 
supremum of the right-hand side of ( |92l ) becomes the integral in (l89l ). ■ 



Now, we prove Theorem [3] following the induction method. 

Proof: [Theorem [3l We first check the Q = 1 case. Note that ti = 1 in this case, thus, (|37] ) implies the 
hypothesis (|38] | holds for Q = 1. 

Next we assume that the hypothesis (1381) holds for Q = j — 1 and consider Q = j. Let 



for m = 1, . . . , J — 1, 



and 



/(c) = -Inc 



/i(c) 



(j 



17 In 122 , 



m=l 



Since < 1 — c, X^^i "^m = induction hypothesis for Q = j — 1 implies 



lim sup Am • P 



i-1 



m'rm 



> C 



,m=l 



< He). 



(93) 
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Note that /(c) and h{c) are monotone decreasing and integrable, and /'(c) is integrable. Then, by Lemma |3j 



lim sup Am P 

Ai,...,A,^oo 



E 

.m=l 



'7~m, 0rn. ^ C 



lim sup JJ Am P 

Ai,...,A,^oo ^^-^ 



i-1 



Tj(/.j + (1 - Tj) 



m=l 



< 



■1 1 



where 



1 r i n 1 

_ C-(l-T,) 



-dz 



(94) 



Since — In z is a convex function , Jensen's inequality implies that 
In 1^ r = — In 



C - ZTj - (1 - Tj - Tm) 



ZTj - C + (1 - Tj) Tm - (1 - C) ^ - ZTj ^ 



1 - C 



+ 



1 -c 



^ ZTj - C + (1 - Tj) Tm - (1 - C) Tj - ZTj ^ 



1 - C 

'^.In. 



\-q C - (1 - Tm) 



Tm 1 - C 

for q< z 1 < < j — 1. 



Hence, 



1-1 i-i 



TVr 



g ^ m=l ^ ^^-J ^-^ ^" 



i-1 



.m=l 



Tra 




z — q 


C - (1 - Tm) 




\-q 



i-1 



(95) 



(96) 



Note that 1/z and [(z — q)j(\ — q)^ ^ are, respectively, monotonically decreasing and increasing in z. Chebyshev 
integral inequaUty [31] implies 

i-i 





Z-q 


Iq Z 


\-q 



dz < 



1 



1 







i-i 


z-q 


dz ■ 1 




l-q 


Jq 



In q 
j 

Finally, by combining ( |94l ). ( |96l ). and ( [97] ). we have the desired result 



(97) 



lim sup Am P 



Ai,...,At — *oo _ -I 
11 1 J m=l 



< 



J-i c- 1 



.m=l 



(98) 



Since both the baseline case (Q = 1) and the inductive step {Q = j) satisfy ( [38] |. we conclude that the hypothesis 
dM]) holds for all Q. ■ 
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